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ABSTRACT 

We outline a method by which the angular radii of giant and main sequence stars located in the Galactic 
Bulge can be measured to a few percent accuracy. The method combines comprehensive ground-based 
photometry of caustic-crossing bulge microlensing events, with a handful of precise (~ 10/ias) astrometric 
measurements of the lensed star during the event, to measure the angular radius of the source, 0*. Dense 
photometric coverage of one caustic crossing yields the crossing time scale At. Less frequent coverage 
of the entire event yields the Einstein timescale t& and the angle </> of source trajectory with respect to 
the caustic. The photometric light curve solution predicts the motion of the source centroid up to an 
orientation on the sky and overall scale. A few precise astrometric measurements therefore yield #Ej the 
angular Einstein ring radius. Then the angular radius of the source is obtained by 0* = 0^{lS.t/t^) sin <f>. 
We argue that the parameters t-E,At,<j), and 6e, and therefore £?*, should all be measurable to a few 
percent accuracy for Galactic bulge giant stars using ground-based photometry from a network of small 
(lm-class) telescopes, combined with astrometric observations with a precision of ~ 10/Ltas to measure 
6*e- We find that a factor of ~ 50 times fewer photons are required to measure £?e to a given precision for 
binary-lens events than single-lens events. Adopting parameters appropriate to the Space Interferometry 
Mission (SIM), we find that ~ 7 minutes of SIM time is required to measure #e to ~ 5% accuracy for 
giant sources in the bulge. For main-sequence sources, #e can be measured to ~ 15% accuracy in ~ 1.4 
hours. Thus, with access to a network of lm-class telescopes, combined with 10 hours of SIM time, it 
should be possible to measure 0» to 5% for ^80 giant stars, or to 15% for ^7 main sequence stars. We 
also discuss methods by which the distances and spectral types of the source stars can be measured. A 
byproduct of such a campaign is a significant sample of precise binary-lens mass measurements. 

Subject headings: stars: fundamental parameters, binaries - gravitational lensing - astrometry 

1. INTRODUCTION 

Although of fundamental importance to stellar astro- 
physics, precise measurements of angular radii are gener- 
ically difficult to acquire routinely and in a model- 
independent way. Classical direct methods of measur- 
ing stellar radii include lunar occultations, interferome- 
try, and eclipsing binaries. Lunar occultation measure- 
ments yield precise angular radii (see Richichi et al. 1999 
and references therein) , but the number of stars to which 
this technique can be applied is limited. The number 
of direct measurements using interferometers has recently 
increased dramatically with advent of, e.g. the Palomar 
Testbed Interferometer (van Belle et al. 1999a, Colavita et 
al. 1999), and the Navy Prototype Optical Interferometer 
(Armstrong et al. 1998, Nordgren et al. 1999), and is likely 
to continue to increase as more technologically advanced 
interferometers come online. Unfortunately, both lunar 
occultation and interferometric angular diameter measure- 
ments have traditionally been primarily limited to nearby, 
evolved stars. Angular radii of main-sequence stars can 
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be determined using detached eclipsing binaries (i.e. Pop- 
per 1998), however the large amount of data (both photo- 
metric and spectroscopic) required to yield accurate radii 
determinations makes this method prohibitive. Thus, of 
the ~ 300 direct, precise angular diameter measurements 
compiled by van Belle (1999), the overwhelming majority, 
~ 85%, are of evolved stars. Finally, it will be difficult to 
acquire a large sample of angular radii determinations of 
stars with metallicity considerably smaller than solar us- 
ing these methods, due to the paucity of metal-poor stars 
in the local neighborhood. 

Here we present a method, based on a suggestion by 
Paczyhski (1998), of measuring angular radii of stars that 
overcomes some of the difficulties inherent in the classical 
methods. This method employs the extraordinary angular 
resolution provided by caustics in gravitational microlens- 
ing events, and as such is yet another in the growing list of 
applications of microlensing to the study of stellar astro- 
physics (see Gould 2001 for a review). The original sug- 
gestion of Paczyhski (1998) was to invert the method of 
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Gould (1994) for measuring the relative source-lens proper 
motion /i ro i in microlcnsing events. If the lens transits the 
source in a microlensing event, precise photometry can be 
used to determine the time it takes for the lens to transit 
one source radius, t* = 0*//z re i, where 9* is angular radius 
of the source. An estimate of 6* , using an empirical color- 
surface brightness relation, together with a measurement 
of the flux of the source, can then be used to estimate /U re i, 
which Gould (1994) argued could be used to constrain the 
location of the lens. However, as Paczyhski (1998) pointed 
out, it is possible to independently measure the angular 
Einstein ring radius of the lens, 
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by making precise astrometric measurements of the cen- 
troid shift of the source during the microlcnsing event us- 
ing, i.e., the Space Interferometry Mission (SIM). 3 Here M 
is the mass of the lens, D is defined by, D = D 0S D i/D} s , 
and D os , D Q \, and D\ s arc the distances between the 
observer-source, observer-lens, and lens-source, respec- 
tively. Since /i rc i = #e/£e, by combining the measure- 
ment of 6*e with the Einstein timescale ts of the event 
determined from the light curve, it is possible to mea- 
sure Osf for the source stars of microlensing events. We 
show that, with reasonable expenditure of resources, it 
should be possible to measure angular radii of a signifi- 
cant sample (~ 80) of giant stars in the bulge to an accu- 
racy of <; 5%, or ~ 7 main-sequence stars to an accuracy 
of <^ 15%. Limb-darkening determinations should also be 
possible for the majority of the sources, and most will be 
relatively metal poor as compared to those for which an- 
gular radii determinations are currently available. 

Although measurements of 0* can be made using single- 
lens events, in §2 we argue that this method is better 
suited to caustic-crossing binary-lens events, which are 
more common, easier to plan for, and considerably less 
resource- intensive than source-crossing single- lens events. 
We describe in some detail the basic method of measur- 
ing 0» for the source stars of caustic-crossing binary-lens 
events in §3, including a discussion of the expected errors 
on the individual parameters that enter into the measure- 
ment. We discuss various subtleties, complications, and 
extensions to the method in §4, and also present an es- 
timate of the number of measurements that might be 
made in this way. Finally, we summarize and conclude in 
§5- 

2. BINARY VERSUS SINGLE LENS EVENTS 



is formally infinite. Large gradients (with respect to source 
position) in the magnification exist near caustics, enabling 
the resolution of the source. Generically, microlenses come 
in two classes: single and binary lenses. Here 'binary lens' 
means a lens systems composed of two masses with an- 
gular separation of order the angular Einstein ring radius 
of the system. Very close and very wide binaries act es- 
sentially as single lenses. Single lenses have a caustic that 
consists of a single point at the position of the lens. In 
these cases, the magnification close to the caustic diverges 
as the inverse of the distance to the caustic. In contrast, 
the caustics of binary lenses are extended, and can cover a 
significant fraction of the Einstein ring. The caustics of bi- 
nary lenses generically consist of two types of singularities, 
folds and cusps. 4 Near a fold, the caustic is well-described 
by generic linear fold singularities, for which the magnifi- 
cation locally diverges inversely as the square-root of the 
distance to the caustic (Schneider & Weiss 1986, Schnei- 
der, Ehlers, & Falco 1992, Gaudi & Petters 2002a). Cusps 
are points where two fold caustics meet, and the magni- 
fication for cusps locally diverges roughly inversely as the 
distance to the cusp point, similar to the magnification 
pattern near the point caustic of a single lens (Schneider 
& Weiss 1992, Schneider, Ehlers, & Falco 1992, Gaudi & 
Petters 2002b). The fact that the magnification near the 
point caustic of a single lens or near a cusp diverges in- 
versely as the distance to the singularity, rather than as 
the square root of the distance to the singularity as with 
folds, means that for a given source size, the 'resolving 
power' of fold caustic crossings is less than single lens or 
cusp crossings. 

Although single-lens events are more well suited to stud- 
ies of stellar atmospheres, they are much less useful for 
measuring the sizes of stars for three main reasons. They 
are generally rarer than fold caustic crossings, their cross- 
ings cannot be predicted in advance, and their centroid 
motion is more difficult to measure. 

Caustic-crossing binaries comprise roughly f cc — 7% of 
all events toward the Galactic bulge (Alcock et al. 2000a; 
Udalski et al. 2000). Of the remaining 1 — f cc events, which 
we will conservatively assume due to single lenses, only a 
fraction #*/#e will exhibit caustic crossings. Therefore, 
the expected ratio of binary-to-single events for which the 
source is resolved (and thus measurement of 9* is possible) 
is 
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The primary requirement to be able to measure #e in 
a microlensing event is that the source should be resolved 
by the gravitational lens. This effectively means that the 
source must cross a caustic in the source plane. Caustics 
are the set of positions in the source plane where the deter- 
minant of the Jacobian of the lens mapping from source to 
lens plane vanishes, and where the magnification therefore 
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4 For binary-lenses, higher-order, beak-to-beak singularities can exist for specific combinations of the binary-lens mass ratio and angular 
separation in units of 9^. However, folds and cusps are the only stable singularities of any lens system (Petters, Levine, & Wambsganss 2000). 
Bcak-to-bcak singularities are unstable in the sense that, for infinitcsimally small variations in the lens parameters, a beak-to-beak singularity 
disintegrates into two cusp-type singularities. Thus the set of parameters where beak-to-beak singularities are expected is formally sparse, and 
practically small. Interestingly, Alcock et al. (2000a) suggested that an observed event may have been due to a source crossing a beak-to-beak 
singularity in a binary-lens, although Alcock et al. (2000b) seem to favor the interpretation that this event is due to a background supernova. 



where i?„ is the physical radius of the source, and we 
adopted D a \ = 6 kpc, D os = 8 kpc, and M = 0.3M Q 
for the scaling relation on the extreme right hand side. 
Although not overwhelming for giant sources, for main se- 
quence sources we expect at least an order of magnitude 
more binary lensing events for which the source is resolved. 
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The ratio of the number of fold-to-cusp crossings is roughly 
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The process of 6** measurement can therefore be subdi- 
vided into three basic steps: 
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where N cusp is the number of cusps, which is cither 6, 8, or 



(1) Measurement of the caustic crossing timescale At 
from a single photometrically well-resolved caustic 
crossing. 

Measurement of the angle <j> and timescale ts from the 
global fit to the binary-lens light curve. 

Measurement of the angular Einstein ring radius #e 
using precise astrometric measurements of the source 
centroid. 
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Y cusp 

10 for a binary lens, depending primarily on d. Thus the 
overwhelming majority of events for which it will be possi- 
ble to measure 9* will be fold caustic crossing binary-lens 
events. 

Binary-lens fold caustic-crossing events also have the 
advantage that the second caustic crossing can be pre- 
dicted in advance. This is because fold caustic crossings 
always come in pairs, and it is typically easy to tell, even 
with sparse sampling, that the first caustic crossing has 
occurred. Then more frequent sampling can be used to 
monitor the rise to the second caustic, in principle en- 
abling the prediction of the time of the second crossing a 
day or more in advance (Jaroszyhski & Mao 2001), and 
allowing the marshaling of the resources necessary to ob- 
tain the dense coverage of the second crossing needed to 
measure the crossing time At of the source (see §3.1). In 
contrast, a single lens caustic crossing can only reliably be 
'predicted' at about the time it begins. 

Finally, it is considerably harder to measure 9^ for 
caustic-crossing single-lens events than binary-lensing 
events. Single lens events have a maximum absolute 
centroid shift relative to the unlcnsed source position of 
whereas binary lensing events can exhibit large 
variations of size ~ #e or more when the source crosses the 
caustic (Han, Chung, & Chang 1999). Therefore consider- 
ably more time will generally be required to determine #e 
to given accuracy for single-lens events than for binary-lens 
events. Since the astrometric measurements essentially re- 
quire the capabilities of SIM (or some similarly precious 
instrument), it is highly desirable to minimize the amount 
of time spent on this step. 

Given the above arguments, we conclude that fold caus- 
tic crossing binary lensing events are the most suitable for 
use in routinely measuring 9*. We will therefore focus on 
this case for the remainder of the discussion, however we 
will briefly revisit single-lens and cusp-crossing events in 
§4.1. 

3. THE METHOD 

Consider a binary-lens event in which the source crosses 
a simple linear caustic. 5 Defining At as the (one-half) 
the time it takes for the source to completely traverse 
the caustic, and (f> as the angle between the source trajec- 
tory and the tangent to the caustic at the crossing point, 
then the time for the lens to cross the angular radius of 
the source is t* = At sin <fi. However, we also have that 
t* = #*//i r ei) where again /i re i is the relative source- lens 
proper motion. Combining these expressions, we have that 
9* = fr c \At sin (j). Using the definition of fj, Te \, we can write 
the angular source radius 9* as the following function of 
observables, 

0, =0E^sin0. (4) 

5 A caustic can generally be approximated as a simple linear fold when the curvature of the caustic is everywhere small on angular scales 
of O(0*), and when the angle of incidence of the source trajectory to the caustic is not small. This approximation will break down when the 
source is large compared to the overall size of the caustic, when the source crosses near a cusp, or when the source 'straddles' the caustic for a 
long time due to a small incidence angle. Although it will still be possible to measure 0-* for such events, the relation between the observables 
and 0* is less straightforward. We discuss such cases in §4.1. 



Each of these steps are illustrated schematically in Figure 
1. 

In the following subsections, we consider each of these 
steps in more detail. We outline the basic requirements for 
the measurement of each of the four parameters (At, </>, t^, 
and 0e), and the expected accuracy with which each can 
be determined assuming reasonable expenditure of observ- 
ing resources. 

3.1. Measuring At 

When the source is interior to a binary-lens caustic, five 
images are created. As the source approaches a fold caus- 
tic, two of these images brighten and merge in a char- 
acteristic way, and eventually disappear when the source 
completely exits the caustic. In contrast to the significant 
brightening of the two images associated with the fold, the 
remaining three images generally vary only slowly over the 
timescale of the crossing. All fold caustics locally have this 
generic behavior, and the magnification A(t) of the source 
near a caustic crossing as a function of time is typically 
well-fitted by the functional form (Albrow et al. 1999b), 



A(t) = 



At 
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+ A cc + uj(t-t cc ), (5) 



where t r is the effective rise time of the caustic, which is re- 
lated to the local derivatives of the lens mapping (Petters 
et al. 2001, Gaudi & Petters 2002), t cc is the time when 
the center of the source crosses the caustic, A cc is the mag- 
nification of all the images unrelated to the fold caustic at 
t = t cc , uj is the slope of the magnification of these images 
as a function of time, and Gq(x) can be expressed in terms 
of complete elliptic integrals of the first and second kind 
(Schneider & Weiss 1987). Note that equation (5) is only 
formally appropriate for a simple linear fold caustic. Thus, 
for a well-sampled fold caustic crossing, At can be deter- 
mined essentially independently of the global geometry of 
the event, and indeed without reference to the photometric 
data away from the crossing itself. In practice, the mag- 
nification is not directly observable, but rather the flux 
F(t) as a function of time. The form for F(t) takes on a 
similar form as equation (5), but with a slightly different 
parameterization (see Albrow et al. 1999b). 
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Fig. 1. — Schematic illustration of the various steps involved in measuring the angular radius of the source star of a caustic-crossing 
binary-microlcnsing event, (a) Infrequent (~ 1 per day) observations of bulge microlcnsing events from a single site reveal a caustic-crossing 
event real-time, enabling a prediction for, and intense photometric monitoring of, the second caustic crossing, (b) Fitting the photometry near 
the second crossing to a generic fold caustic model yields the caustic crossing timescale At. (c,d) Fitting to the entire photometric datasct 
to a binary-lens model yields the Einstein timescale *e of the event, as well as the angle cf> of the trajectory with respect to the caustic, (e) 
The photometric solution predicts the motion of the centroid of the source up to an unknown scale and orientation on the sky. A few, precise 
astrometric measurements can then be used to determine 6*e, the angular Einstein ring radius. The cross shows errorbars of lO^tas. The 
angular source size is then given by 6* = d-^{At/t-g,) sin0. 
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Equation (5) assumes a uniform source. This will likely 
be a poor approximation in optical bands, and assuming 
uniform source in the presence of limb-darkening may re- 
sult in a systematic underestimate of At, and therefore 
since, the effect of limb-darkening can partially com- 
pensated for by a smaller (dimcnsionless) source size, at 
least for poorly-sampled light curves. However, for well- 
sampled caustic crossings and <~ 1% photometric accuracy, 
it should generally be possible to accurately measure both 
the source size and limb-darkening coefficient (s) (Rhie & 
Bennett 1999). Such data quality is readily achievable; 
indeed, independent limb-darkening and (dimensionless) 
source size measurements have been made for the source 
stars of at least five microlensing events (Albrow ct al. 
1999a, Afonso et al. 2000, Albrow et al. 2000, Albrow et al. 
2001, An et al. 2002). A generalized form of equation (5) 
that includes simple limb-darkening can be found in Al- 
brow et al. (1999b) and Afonso ct al. (2001). 

By fitting equation (5) (or a generalized form of it) to 
the light curve near a well-sampled fold caustic crossing, 
one can derive the parameters t r , t cc , A cc , uj, and At. The 
parameters t cc , A cc , and t r can subsequently be used to 
constrain the global solution to the entire light curve (see 
Albrow et al. 1999b). However, of primary interest here 
is the parameter At, whose value is essentially indepen- 
dent of the global solution. See Figure 1(b). This means 
that At and the parameters determined from the global 
solution, t~E and <fi, will be essentially uncorrelated. 

In order to be able to determine At from the caustic 
crossing data alone, the caustic crossing must be well- 
sampled, so that the parameters in equation (5) be well- 
constrained. Practically, this requires forewarning of the 
caustic crossing. Fortunately, this is generally possible 
with only photometry available from the collaborations 
which survey the Galactic bulge and find the alert the 
microlensing events real-time (EROS, Afonso et al. 2001; 
MOA, Bond et al. 2001; OGLE, Udalski et al. 2000^ 
Woziak ct al. 2001), although improved predictions would 
be possible with continuous photometry (Jaroszynski & 
Mao 2001). Thus no additional resources need to be in- 
vested to predict caustic crossings. However, as we discuss 
in §3.2, some additional sampling of the overall light curve 
may be needed to constrain t^ and <p an d determine a 
unique global solution. 

The accuracy with which At can be determined for a 
given caustic crossing will depend not only on the intrinsic 
parameters of the caustic crossing, but also on the sam- 
pling rate and photometric accuracy near the caustic cross- 
ing, which in turn will depend on weather, blending, etc., 
which tend to vary in a stochastic manner. Therefore it is 
not very useful to attempt to quantify the expected errors 
on At for an idealized observing setup. However we can ob- 
tain an order-of-magnitude estimate for a^t by examining 
measurements of At from published analyses of observed 
caustic-crossing events. We discuss these determinations 
more thoroughly in §3.4. Typically, well-covered caustic 
crossings yield fractional errors of a At/ At <~ 1%. 

3.2. Measuring tE and <p 

The parameters tE and 4> are determined from the global 
solution to the overall geometry of the binary-lens light 
curve. The relationship between the salient features of 
an observed light curve, and the canonical parameters of a 



binary- lcnsing event, are generally not obvious or straight- 
forward. This fact generally makes fitting an observed 
light curve, and thus inferring the parameters tE and <f>, 
quite difficult (see Albrow et al. 1999b for a thorough dis- 
cussion) . Although many methods have been proposed to 
overcome these difficulties, the lack of obvious correspon- 
dence between these parameters of interest and the light 
curve morphology generally implies that it is difficult to 
make general statements about the kinds of observations 
that are needed to reliably measure tE and <f>. However, 
we can make some generic comments. For caustic-crossing 
binary lenses, one potential observable is the time St cc be- 
tween caustic crossings. Between caustic crossings, the 
magnification is typically considerably larger than outside 
the caustic, and therefore even with sparse (~ 1 day) sam- 
pling, it should be possible to determine to reasonable 
precision the time of the first caustic crossing, provided 
that data exist before the first crossing. The requisite 
photometry will generally be acquired by the survey col- 
laboration^). Once the binary- lens event is alerted, more 
frequent photometry can be acquired by follow-up collab- 
orations (PLANET, Albrow et al. 1998; MPS, Rhie et al. 
2000), thus mapping the "U" -shaped curve between the 
caustic crossings. This shape, combined with information 
from cusp-approaches (or lack thereof) just outside the 
caustic, provides information about the shape of the caus- 
tic and the trajectory of the source through it. From this, 
the angle <p of the trajectory with respect to the caustic 
can be derived, and also the distance d cc between the caus- 
tic crossings in units of 6>e- Then, the timescale is given 
by ts = t cc /d cc . 

Although the above analysis is highly trivialized, it does 
suggest that the following steps should be taken to en- 
sure an accurate measurement of ts and <f>. First, it is 
important to constrain the time of both caustic crossings 
to reasonable precision. This means that the survey col- 
laborations should sample on timescales no less than a few 
days. Second, followup photometry should be initiated rel- 
atively soon after the first crossing, to measure the shape 
of the intra-caustic light curve reasonably well. This is 
also necessary in order to predict the second caustic cross- 
ing (Jaroszynski & Mao 2001). Furthermore, the follow-up 
photometry should continue past the second caustic cross- 
ing, to detect cusp-approaches or cusp-crossings (or the 
lack thereof). 

Due to the difficulties inherent in fitting binary-lens light 
curves, it would be extremely difficult to attempt to pre- 
dict the expected errors on ts and 4> for a hypothetical 
observing scenario. Furthermore, due to the complicated 
relation between observables and parameters implies that 
these errors are likely to depend strongly on the geometry 
of the event and light curve coverage, and therefore such 
'predictions' would not be very useful. However, we would 
like to have an order-of-magnitude estimate for the ex- 
pected errors given reasonable light curve coverage. From 
determinations of these parameters in published caustic- 
crossing events, we can expect fractional accuracies of a 
few percent, provided that the light curve is well-covered, 
in the sense outlined above. However, if the light curve 
coverage is incomplete, then errors of ;> 20% are expected. 
See §3.4. 

3.3. Measuring 6*e 
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Event Name 


f At/ At 


0-0/0 (0) 




C* E /tE 


I 


Rcf. 


Comments 


MACHO 98-SMC-l 


0.9% 


32% (36°9) 




23% 


22.1 


1,2 


Fold crossing 


OGLE-1999-BUL-23 


0.5% 


0.2% (56°1) 




1% 


18.1 


3 


Fold crossing 


MACHO 95-BLG-30 






0.1% 




13.4 


4 


Single lens 


MACHO 97-BLG-28 






0.3% 




15.6 


5 


Cusp crossing 


MACHO 97-BLG-41 






5% 




16.8 


6 


Rotating binary 


EROS 2000-BLG-5 






0.8% 




16.6 


7 


Parallax effects 



Table 1 Fractional Errors on Parameters from Observed Events. 



References.— (1) Albrow ct al. (1999b), (2) Afonso ct al. (2000), (3) Albrow et al. (2001a), (4) Alcock et al.(1997), (5) Albrow et al.(1999a), 
(6) Albrow ct al. (2000), (7) An ct al.(2002) 



A global solution to the entire light curve effectively re- 
quires the specification of the vector position of the source 
u(t) as a function of time in units of 9e, and the topology 
of the lens, i.e. the mass ratio q and projected separation 
d in units of 9e- These parameters yield not only the total 
magnification A of all the images as a function of time, but 
also the individual image positions <f i and magnifications 
Ai as a function of time, Therefore, it is also possible to 
predict the centroid ip cl of all the individual microimages, 



The accuracy with which 9e can determined will depend 
on the geometry of the event, the time of the astrometric 
measurements, and the time span between the measure- 
ments. Furthermore, as we discuss, there are, in reality, 
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Note that ip ci is the centroid shift with respect to the lens 
position. Is it customary to consider the centroid shift with 
respect to the unlensed source position, 5ip cl = (p cl — u. 
The unlensed source position it, which is comprised of the 
parallax and proper motion of the source in some astro- 
metric frame, can be determined via measurements of the 
unlensed motion of the source. Since the astrometric ef- 
fects falls off very slowly (as these must be obtained 
many tg after the event is over. In fact, such measurements 
are not strictly needed in order to measure 9e (but may 
be desirable for other reasons, see §4.3). Rather, one can 
simultaneously fit astrometric measurements during the 
course of the event for the relative position and proper mo- 
tion of the source (in order to establish a local astrometric 
reference frame), and the offset induced by microlensing. 

The centroid shift (p cl is in units of 0e, and its compo- 
nents are oriented with respect to the projected binary- 
lens axis, whose orientation a on the sky is unknown. The 
observable centroid is 
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Fractional Error o g /9 E 

FIG. 2.- The cumulative distribution of expected errors in Be in a 
Monte Carlo simulation, assuming a contribution of blended light of 
0% (no blending, solid), 1% (dotted), 10% (long dashed), and 90% 
(dot-dashed). See text. 

additional parameters that must be determined from the 
astrometric data. We therefore perform a Monte Carlo 
simulation to estimate the accuracy with which 6e can 
be recovered. We explore whether the various parameters 
can be measured independently, or are degenerate with 
the measurement of 9e, and study the range of fractional 
uncertainties in measuring 9e for a full ensemble of binary 
lenses. 

Our simulation consists of the following elements: a 
Galactic model and model of the ensemble of lenses to 
generate an ensemble of microlensing events, an observa- 
tional strategy, and a fit of these observations to a set of 

technically, this is only true for an astrometric observer that is spatially coincident with the photometric observer. This will not be true 
for astrometric observations with SIM which will be in an Earth-trailing orbit. We discuss this complication in §4.2.2. 



Alternatively, one can combine a and 6>e and simply con- 
sider the vector Oe- Thus the global solution to the photo- 
metric light curve yields not only and <f>, but also yields 
a prediction for the astrometric curve, up to an unknown 
orientation and scale 6*e- 6 Therefore, by making a series of 
astrometric measurements 9 c i (t) at several different times 
during the course of the event, one can determine 9e and 
a via equation (7), and using the predictions for y c ;(t) 
from the photometric solution. 
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parameters including Oe- 

Our ensemble of lenses is identical to that of Graff & 
Gould (2002). Briefly, we draw sources and lenses from 
sclf-lensing isothermal sphere 8 kpc from the observer, and 
with two-dimensional velocity dispersion of 220km s _1 . 
Both masses in the binary lens are chosen from the rem- 
nant mass function of Gould (2000). We pick a flat distri- 
bution in log(d), the logarithm of the dimensionless binary 
separation. The path of the source through the lens ge- 
ometry is chosen randomly with a uniform distribution of 
angular impact parameter We, and we only consider paths 
that cross caustics. 7 

As discussed in Graff & Gould (2002), in this ensem- 
ble of events there are many more events with a short 
time between caustic crossings, the caustic interior time 
ti n t, than is observationally detected by the MACHO and 
OGLE observing groups. This led these authors to sug- 
gest that most events with a short caustic interior time 
are not detected as caustic crossing binaries, and to define 
a caustic crossing detection efficiency £i,(fi n t) akin to the 
standard single lens detection efficiency £ s (tE)- 

The observational strategy is relatively unimportant, as 
long as there are astrometric measurement on either side 
of a caustic. Although it is conceivable that observations 
might be scheduled at particularly favorable times, such 
as times of maximum magnification or maximum displace- 
ment of the image centroid, it is likely that the telescope 
measuring the astrometric displacement will be oversub- 
scribed. A simpler strategy would be to schedule periodic 
observations in advance. We have assumed that observa- 
tions will be made every four days, with a 24 hour delay 
after the event is recognized as a caustic crossing binary, 
i.e., after the first caustic crossing. That is, the first obser- 
vation comes 1 — 5 days after the first caustic crossing. We 
assume observations are made for a total of 36 days, which 
corresponds to approximately 2tE for the median event 
timescale. Practically, observations should continue until 
after the second caustic crossing, and the total number of 
observations should be at least as large as the number of 
parameters to be constrained. 

As we discuss in §3.4 (see also Table 1), well covered bi- 
nary lenses can be fit photometrically with small errors on 
the parameters. Thus, we have assumed that all the pa- 
rameters which can be fit from a single photometric tele- 
scope are determined. 

Given a binary microlensing event from our ensemble, 
we use our observational strategy to create a series of pho- 
tometric and astrometric measurements, respectively A(t) 
and c i (t) which we can combine into a single list of mea- 
surements Mfc each with uncertainty Ok ■ Using the Fisher 
matrix technique, e.g., Gould & Welch (f996), we deter- 
mine the co variance matrix of the errors 
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and that the total telescope time, aperture, efficiency, fil- 
ter width, and source brightness are such that a total 
of N — 60, 000 photons would be detected from an un- 
magnified source for a total photometric signal to noise of 
250. 8 We assume that the fractional photometric accuracy 
is simply iV -1 / 2 , and that the astrometric uncertainty is 
ag = N~ 1 / 2 0f, where Of is the width of the point spread 
function, or in the case of an interferometer, the fringe sep- 
aration. Here we have assumed Of — 2.5mas. It is trivial 
to scale our results to brighter sources or larger telescopes: 
the fractional uncertainty in 6e is simply proportional to 

AT-l/2. 

We always fit for the four parameters required to estab- 
lish a local astrometric frame, and 0e and a, the size scale 
and orientation of the microlensing excursion. We also as- 
sumed that the satellite which measures the astrometric 
motion is 0.2 AU from the (ground based) photometric 
measurements which fix the lens parameters. Thus, we 
can simultaneously fit for f^, the projected Einstein ring 
radius, in the manner of Graff & Gould (2002). In ad- 
dition to this basic fit, we also considered blending from 
luminous lenses and binary sources, which requires several 
additional parameters. We discuss parallax and blending 
in §§4.2.2 and 4.2.3, respectively. 

Our basic results are summarized in the leftmost curve 
in Figure 2. We see that in the absence of blending, 0e can 
be determined with < 10% uncertainty in 85% of events. 
The median error is cto^/Oe — 2.4%. This is comparable to 
the uncertainty found by Gould & Salim (1999) for single 
star events, but with 50 times as many photons as we have 
assumed here. Thus, it is much easier to measure 0e for 
caustic crossing binary lenses than for single lenses. 

3.4. Expected Fractional Errors on 0* 

From the expression for 0* (Eq. 4), and assuming the 
errors in 0e, <e, At, and <j> are small and uncorrelated, the 
fractional error in 0* is given by, 
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Here the ai are the various parameters being fit. The error 
in parameter a» is simply a ai — ^Jcu- 

We assumed that the photometric uncertainty, Ufe, of 
the interferometric telescope, is photon-noise-dominated, 

7 In contrast to the usual technique for single lenses in which events are chosen from a uniform distribution in b, but arc weighted towards 
large $e events by multiplying the mass function by M 1 / 2 . 

8 Note that a telescope with an overall efficiency of ~ 30% and diameter Aj- collects ~ 100(Ay/lm) 2 photons per second at I = 18. 



(9) 

where ag E1 o~t E , o~At, and are the uncertainties in 0e, 
tE, At and </>, respectively, and is in radians. 

To date, there exist in the published literature 10 mi- 
crolensing events for which the source star was well- 
resolved. Unfortunately, for four of the events, those pre- 
sented in Alcock et al. (2000a), no estimate of the errors of 
the derived parameters is given. We therefore cannot use 
these events to explore the expected magnitudes of a tE , 
a A t, and a^. 

For the six events for which errors on the relevant fit 
parameters were given or derivable, only two of them are 
generic binary-lens fold caustic crossing events. The er- 
rors <7t E , (j At; and (70 for these two events, MACHO 98- 
SMC-I and OGLE-1999-BUL-23, are presented in Table 
1. For MACHO 98-SMC-l, these errors have been deter- 
mined from the ensemble of solutions presented in Albrow 
et al. (1999b), which were fits to the PLANET collabora- 
tion photometry, which only covered the last half of the 
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event, whereas the analysis of the combined photometry of 
the EROS, MACHO, MPS, OGLE, and PLANET collab- 
orations (Afonso et al. 2000) yields considerably smaller 
errors. We concentrate on the results of Albrow et al. 
(1999b) here in order to demonstrate the kinds of errors 
that result from incomplete light curve coverage. Note 
that a At ~ 1%, which is not surprising, since the accu- 
racy with which At can be determined depends almost 
exclusively on the photometric coverage near the caustic 
crossing. However, the errors on the global parameters £e 
and (j> are quite large, a tE ~ 20% and ~ 30%. This is 
due to the fact that the PLANET photometry only cov- 
ered the latter half of the event, and contained no data 
prior to or during the first caustic crossing. Therefore the 
global geometry was quite poorly constrained with their 
data alone. Such incomplete coverage would clearly jeop- 
ardize a precise measurement of 0*. This is in contrast to 
OGLE-1999-BUL-23, for which both PLANET and OGLE 
obtained data before the first caustic crossing. In this case, 
it was possible to measure At, t~E, and <j> to better than 
1% (Albrow et al. 2001a). 

The remaining four include a single lens event (MA- 
CHO 95-BLG-30; Alcock et al. 1997), a binary-lens event 
in which the source crossed a cusp (MACHO 97-BLG-28; 
Albrow et al. 1999a), a binary- lens event in which the ro- 
tation of the binary was detected (MACHO 97-BLG-41; 
Albrow et al. 2000), and a binary-lens event for which 
both rotation and parallax effects were detected (EROS 
BLG-2000-5; An et al. 2002). In general, the previous dis- 
cussion is not directly applicable to these events, as the 
information in the light curves is not easily decomposed 
into the parameters At, (j> and £e- Nevertheless, in all four 
cases the dimensionless source size = 0*/0e was deter- 
mined. In three of the cases, <5/9*/p* <; 1%. In one case, 
MACHO 97-BLG-41, <5p*/p* - 5% , due primarily to the 
fact that there exists only a handful of data points in which 
the source was resolved. Thus, although it is difficult to 
draw any general conclusions from these unique events, it 
is does seem likely that errors of <J 1% are achievable for 
most types of events in which the source is resolved. 

Also presented in Table 1 are the determined values of 
the /-magnitude of the source, for all six events. At I <~ 22, 
the source star for MACHO 98-SMC-l would be too faint 
to target with most upcoming interferometers, including 
SIM. The other events are primarily bulge clump giants 
(I ~ 15), for which accuracies of <7e E /#E <J 5% should be 
achievable with a reasonable amount of exposure time with 
upcoming interferometers, and in particular with <; 8 min- 
utes of SIM time (see §4.5 for a estimate of the required 
exposure times for SIM). In these cases, the expected error 
on 0* is typically dominated by ag E , and therefore we can 
expect a e J9* ~ crg E /0 E ~ 5%. The source of OGLE-1999- 
BUL-23 is G/K subgiant (i?» ~ 3i? Q ), and thus dimmer. 
The time required to achieve an accuracy of o~$ e /0e ~ 5% 
will be larger by a factor of ~ 16, or ~ 2 hr for SIM. 

4. DISCUSSION 

Our goal in §3 was to capture the essence of the method 
of measuring 0* , and the discussions were therefore some- 
what oversimplified, and glossed over several important 
points. In particular, we concentrated on fold caustic 
crossing binary-lens events toward the bulge, whereas mea- 



surements of should be possible in other, rarer, types of 
events, such as cusp-crossings and single-lens events, and 
possibly events toward the Magellanic Clouds. We also 
ignored various higher-order effects which could, in prin- 
ciple, complicate the measurements. We therefore briefly 
discuss some of these complications and extensions. We 
also discuss the prospects for measuring the spectral type 
of the source, and also its distance, in order to convert 
from angular radius 0» to physical radius Finally, we 
present an example observing campaign aimed at measur- 
ing angular radii for a significant sample of sources, out- 
lining the resources required, and estimating the number 
of 9* measurements that might be made per year for such 
a campaign. 

4.1. Single-lens, Cusp- crossing, and Magellanic Cloud 

Events 

Although we have focussed on fold caustic-crossing 
binary-lens events toward the bulge, it is important to 
emphasize that 9* can, in principle, be measured for 
other types of caustic-crossing events such as cusp-crossing 
events, single-lens events, and all types of caustic-crossing 
events toward the Magellanic clouds. Indeed, in §3.4 we 
discussed examples in the literature of a single-lens and 
two cusp-crossing events for which a ~ 5% measurement 
of would have been feasible. 

In general, isolated cusp crossings, such as in MACHO 
98-BLG-28, cannot be predicted in advance, and thus 
planning for such events is difficult, if not impossible. 
However, one will still have some advance warning of those 
cusp events which occur just after or in place of second fold 
caustic crossings. For such events, sufficient photometric 
coverage of the crossing should routinely be possible. In 
all cases, it is more difficult to disentangle the information 
arising from the cusp itself with the information from the 
global light curve. This generally implies that the analysis 
of these light curves will be more complicated, however this 
does not necessarily preclude an accurate measurement of 
9^. 

Single-lens events are less desirable simply because they 
require a factor of ~ 50 times more astrometric observ- 
ing time to achieve the same fractional accuracy in #e as 
binary-lens events. Since the astrometric observations are 
likely to be the most limited resource, this makes single 
lens events considerably less attractive. 

If it were possible to measure angular radii of stars in the 
Magellanic clouds (MCs), this would be quite interesting, 
due to the metal-poor nature of the stars. Unfortunately, 
there are several major hindrances to measuring 9* for a 
substantial number of stars in the MCs. First, the event 
rates toward both the MCs are small, and a large number 
of stars must be monitored just to detect a few events per 
year. Therefore, the number of caustic-crossing events is 
quite low. To date, there have been only two caustic cross- 
ing events toward the MCs: MACHO 98-SMC-l, which 
we discussed in §3.4, and MACHO LMC-9. These events 
have source magnitudes of Vs = 22.4 (Afonso et al. 2000) 
and Vs — 21 A (Alcock et al. 2000a), respectively, which 
brings up a second difficulty: SIM cannot follow source 
stars fainter than V ~ 20, so these two events could not 
have been used to measure the angular radii of their source 
stars. In fact, even if the entire LMC were monitored for 
microlensing, only ~ 1 event per year would have Vs ^ 20, 
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and this event would be from an evolved star. The prob- 
ability of a caustic-crossing event (either binary or single 
lens) is smaller by at least an order of magnitude. The 
paucity of events and faintness of the source stars might 
be circumvented if sufficiently rapid Target of Opportunity 
times are available. In this case, it might be possible to use 
intrinsically fainter source stars, for which caustic-crossing 
event will be more common, and measure the astromet- 
ric displacement during the brief period of time when the 
source is highly magnified as it crosses the caustic. The 
maximum magnification of a source of dimensionless size 

— 1/2 

p* crossing a fold caustic is A max <~ p* ' . For main- 
sequence sources, A m&K ;> 30, or more than three magni- 
tudes, and thus sources with Vg <J 23 can briefly be bright- 
ened to SIM detectability. For example, the source star of 
MACHO 98-SMC-l was brighter than V ~ 18 for about 
7 hours during the second caustic crossing. Finally, even 
if the source does attain a sufficient brightness to be mea- 
surable by SIM, it remains to be seen whether the centroid 
varies sufficiently during this time to provide an accurate 
measurement of 9e- This is especially difficult in light of 
the fact that typical value of 9e for self-lensing events to- 
ward the MCs are only an order of magnitude larger than 
SIM's accuracy (Paczyhski 1998, Gould & Salim 1999). In 
summary, it appears that it will be quite difficult to mea- 
sure angular radii of stars in the MCs using this method, 
especially if the majority of the events seen toward these 
targets are due to self-lensing (Sahu 1994). 

4.2. Complications to 9e Measurement 

The method we have presented here is only interesting 
if it can feasibly be used to make precise measurements 
for a large number of sources with reasonable expenditure 
of resources. Since the requisite astromctric instruments 
are likely to be the most limited resource, it is crucial that 
accurate and unambiguous determinations of 9e be generi- 
cally possible using a few astrometric measurements, when 
combined with the photometric light curve solution. We 
have explained how a complete photometric solution gen- 
erally leads to a prediction for the astrometric centroid 
shift up to an unknown scale 9e and orientation a on the 
sky. However, this is true only under a number of simplify- 
ing assumptions, including uniform motion of the observer, 
source, and lens, dark lenses, isolated sources, and unique 
global solutions. If one or more of these assumptions are 
violated, then the prediction for shape of the astromet- 
ric curve may not be unique, and thus the measurement 
of #e may be compromised. We therefore discuss each of 
these complications and under what conditions they may 
be important. 

4.2.1. Binary Lens Degeneracies 

Binary lenses are characterized by two quantities: q, 
the mass ratio, and d, the instantaneous projected sep- 
aration in units of 9e- It has been demonstrated both 
theoretically (Dominik 1999a) and observationally (Afonso 
et al. 2000, Albrow et al. 2002) that certain limiting cases 
of binary lenses can exhibit extremely similar observable 
properties. In particular, Dominik (1999) showed that 
the binary-lens equation can be approximated by an sin- 
gle lens with external shear, or Chang-Refsdal (CR) lens 
(Chang & Refsdal 1979,1984), near the individual masses 



for widely-separated binaries (d 3> 1), and near the sec- 
ondary (least massive) lens when d <C 1. Furthermore, 
near the center-of-mass of a close binary, the lens equation 
is well-approximated by a quadrupolc lens, and both the 
quadrupole lens and CR-lens can exhibit extremely similar 
magnifications when the quadrupole moment is equated to 
the shear (Albrow et al. 2002). Thus there can exist multi- 
ple degenerate solutions to an observed photometric light 
curve, even with extremely accurate photometry. However 
the astrometric behavior of these degenerate solutions is 
very similar both in the shape and overall scale of the as- 
trometric curves, at least for the close/wide degeneracy 
(Gould & Han 2000). Therefore this degeneracy should 
not affect the determination of 9e using the prediction 
from the light curve. It is likely that the other intrin- 
sic denegeracies will also not affect the determination of 
Be, since the degeneracy arises from the lens equation it- 
self, and thus affects both the photometric and astrometric 
curve in the same manner. 

Note that it is important that the normalization of 9e be 
consistent for the two degenerate solutions. For example, 
consider the case of the close/wide degeneracy in MACHO 
98-SMC-l (Afonso et al 2000, Gould & Han 2000). If one 
normalizes to the total mass of the binary, the close solu- 
tion implies a value of 9e, c = 76/ias, whereas the wide so- 
lution has 9e, w = 170yuas. Since the astrometric curves are 
essentially identical (both in shape and scale), one might 
therefore suspect the inversion of this process would yield 
two equally likely values of that differed by a factor 
of <~ 2. Of course, this 'ambiguity' is wholly artificial, 
and arises because the value of 9e, w for the wide binary is 
normalized to the entire mass of the binary, whereas the 
lensing effects are basically caused by the least massive 
lens, since d — 3.25. Normalizing to the mass of the sin- 
gle lens, 9' Ew = 9e. w {1 + q~ 1 )~ 1 / 2 , where here q = 0.24, 
and thus 9 E — 75/ias, essentially identical to the close- 
binary solution. Note that as d approaches unity, the iden- 
tification of the 'proper' 9e normalization becomes more 
nebulous, since the lenses can no longer be considered in- 
dependent. However, the degeneracies also become less 
severe as d —> 1 . 

Dominik (1999b) has also shown that poorly-sampled 
binary lens light curve can also yield distinct degenerate 
solutions. Note that these solutions are 'accidental' in 
the sense that they do not arise from degeneracies in the 
lens equation itself. Thus Han et al. (1999) found that 
such degenerate photometric light curves yield astromet- 
ric curves which are widely different. Such degeneracies 
would prohibit the measurement of 9e using a few astro- 
metric measurements. Therefore well-sampled photomet- 
ric light curves are essential for reliable measurements of 
0*. 

4.2.2. Parallax 

If the two observers are displaced by a significant frac- 
tion of Pe = D9e, the angular Einstein ring radius pro- 
jected onto the observer plane, then the source position 
u as seen by the two observers will be significantly differ- 
ent. Since SIM will be in an Earth-trailing orbit, it will 
drift away from the Earth at a rate of ~ 0.1 AU yr _1 . 
Thus after 2.5 years (half way through the SIM mission), 
it will be displaced from the Earth by t <~ 0.25 AU, which 
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corresponds to a displacement in the Einstein ring of, 



\6u\ = — I sin-yl, (10) 

where 7 is the angle between the the line of sight and the 
Earth-SIM vector. For typical bulge parameters, 



f E = 7.6 AU 



M 
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and therefore \Su\ ~ 3%. This implies that both the mag- 
nification A and the centroid shift 56 c i will be significantly 
different (at a given time) as seen from the Earth and SIM. 
Since the value of te is not known a priori, \5u\ cannot be 
predicted from the ground-based photometry alone, and 
must be estimated from the astrometric data itself. For- 
tunately, SIM will likely have excellent photometric capa- 
bilities (see Gould & Salim 1999), and thus the relative 
magnifications between the light curves from the ground 
and SIM will provide additional constraints. Indeed, in 
our Monte Carlo simulations, we assumed that the astro- 
metric observer was displaced by 0.2 AU from the Earth, 
and simultaneously fit for both #e and te (among other 
parameters), and found that # E could still be constrained 
quite accurately. This is because the information about 
fs comes primarily from the photometry, while the in- 
formation about 6e comes primarily from the astrometry. 
Therefore, the two parameters are not degenerate. Note 
that a 'byproduct' of these measurements is a determina- 
tion of the total mass of the binary lens (Gould & Salim 
1999; Han & Kim 2000; Graff & Gould 2002) 

There are two additional parallax effects. One is due 
to the motion of the Earth (or SIM) around the sun, and 
will become significant on timescales that are a substan- 
tial fraction of a year, which corresponds to many £ E for 
typical bulge events. There is also a second order effect 
that arises from the difference in projected velocities be- 
tween the Earth and SIM. This effect is cx v®/v, where 
v = f^/t-E is the transverse velocity of the lens projected 
on the observer plane, and is ~ 800 km s _1 for typical 
bulge self-lensing events. Since the velocities and posi- 
tions of the Earth and SIM will be known, both of these 
effects can easily be included in the fit for the microlensing 
parallax, and so do not present any additional difficulties. 

4.2.3. Luminous Lenses and Binary Sources 

With its planned 10 meter baseline, SIM will have a 
resolution of ~ lOmas, sufficient to resolve the majority of 
unassociated nearby stars that are blended with the source 
in ground-based photometry (Han & Kim 1999). Since 
the photometric blending is well-constrained in binary-lens 
events, unambiguous prediction of the unblended astro- 
metric behavior of the source is possible. Thus, blending 
will typically not affect the measurement of #e- However, 
luminous lenses and companions to the source star with 
separations <; 10 mas will not be automatically resolved 
by SIM (Jeong, Han & Park 1999). Dalai & Griest (2001) 
have shown that, using two pointings, this limit may be 
lowered to ~ 3mas, however it is essentially impossible to 
resolve multiple sources with separations below this limit 
(e.g. binary source companions). In these cases, all that 



will be measured is the total centroid of all the sources in 
the resolution element. 

The centroid in the presence of luminous lenses, <p cLb , 
is related to the centroid in the absence of blending, ip cl , 

by 



Vcl.b 



<Pcl 



ft 



-<fb 



1 



(12) 



where ft = XI- F^j/Fp is the sum of the flux of all un- 
lensed sources (blends) relative to F , the baseline flux of 
the lensed source, A is the magnification of the source, 
and (p b is the centroid of the blends relative to the ori- 
gin of the lens. From equation (12), it is clear blending 
is more complicated in astrometric microlensing than in 
photometric microlensing: whereas photometric blending 
can be described by one parameter, the blend fraction f t , 
astrometric blending requires two additional parameters, 
the centroid of light of the blend ip b . A special case of 
astrometric blending is bright lens blending. In the single 
lens case, this eliminates the blend location parameters, 
the location of the centroid of light is the moving lens (i.e. 
tp b = 0). In bright binary lens blending, only one parame- 
ter is eliminated, the centroid of light is somewhere on the 
lens axis between the two stars in the lens. However, it will 
generally not be known a priori which case one is dealing 
with, and therefore (p b must be included as a parameter 
in the astrometric fit. 

Blending is problematic because it effectively 'dilutes' 
the astrometric shift between two points in the lightcurve, 
which is qualitatively similar to the effect of changing E . 
If the event is not well covered, these two effects can 
be quite degenerate. In order to determine how degen- 
erate blending is with the E , we have included in our 
Monte Carlo simulations a fixed amount of blended light 
of ft = 1%, 10%, and 90%. We assume that ft is known 
(i.e. from photometry), but <p b is not. The results are 
shown in Figure 2. We find that, if the blending fraction 
is close to 1, then the two effects are nearly degenerate, 
and our fractional uncertainty in E increases by two or- 
ders of magnitude to of order unity. However, for f t <, 1%, 
the median error increases by less than a factor of two. In 
most cases, the blending will be known to be small from 
the photometric data. In these cases, the fractional un- 
certainty in 8e will not be seriously degraded. The few 
events with known large blending can be easily jettisoned 
from the sample. 

4.2.4. Lens Rotation 

The photometric effects of lens rotation in binary mi- 
crolensing events has been explored theoretically by Do- 
minik (1998) and has been detected in event MACHO 
1997-BUL-41 (Albrow et al. 2000). The astrometric effects 
of rotating binary-lenses have not been explored, and it is 
therefore difficult to draw any general conclusions as to 
the importance of this effect. However, to the extent that 
it is detectable in the photometric light curve, lens rota- 
tion poses no difficulties, as its astrometric effect should 
be predictable from the global solution. Effects that are 
photometrically undetectable but astrometrically signifi- 
cant are potentially problematic. 

The amount that a binary-lens rotates during £e is given 
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by, 



ip ~ 4.5°cf 



-3/2 



M 



0.3M Q 



1/4 



1.5kpc 



-1/4 



V 150km s 

(13) 

where D' = D Q \D\ s j D os , and assuming circular, face-on 
orbits. Because the caustic cross-section is maximized for 
binaries with separations of order 0e, the majority of de- 
tected caustic-crossings will have d ~ 1 (Baltz & Gondolo 
2001). Therefore, for typical events, the effects of binary 
rotation should be small if astrometric observations are 
closely spaced with respect to the event timescale £e- This 
is generally also advantageous for the accurate recovery of 
9 E (See §3.3). 

4.3. Measuring D os 

In this paper, we have emphasized the measurement of 
the angular radius rather than the physical radius, 
However, it may also be interesting to measure i?* for some 
events. In order to do this, the distance to the source star 
must be measured independently. Fortunately, the astro- 
metric accuracy needed to measure 6** is generally suffi- 
cient to measure the parallax ir s of the source stars, 



resolution element can be achieved with exposure times of 
tens of minutes for low-resolution spectra (Lennon et al. 
1996), or a couple of hours for high- resolution spectra 
(Minniti et al. 1998), using 8m or lOm-class telescopes. 
Although low-resolution spectra are sufficient for accurate 
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In order to measure i?* to a similar accuracy as (~ 5%), 
ir s must be measured to somewhat better accuracy, which 
implies an astrometric error of erg <J 5/ias. For SIM and 
an / = 18 source, this is achievable with ;> 7 hours of in- 
tegration, which is considerably more time than is needed 
for the 6** measurement alone. However, it is important to 
note that these measurements can be made after the mi- 
crolensing event is over. Therefore, it should be possible 
to employ ground-based interferometers for the measure- 
ment of 7r s , rather than spend precious SIM resources, at 
least for brighter sources. 

4.4. Typing the Source Star 

A measurement of is essentially useless if the spec- 
tral type and luminosity class of the source is not known. 
The source stars of microlensing events can be typed in 
two ways. The first is to simply measure the color and 
apparent magnitude of the source. This information is 
generally acquired automatically from the fit to the photo- 
metric data of the microlensing event. 9 By positioning the 
source star on a color-magnitude diagram of other stars in 
the field, one can generally type the source to reasonable 
accuracy. The primary pitfalls of this method are differen- 
tial reddening and projection effects (i.e. the source may 
in the foreground or background of the bulk of the stars 
in the field). 

A more robust way of typing the source star is to acquire 
spectra. This is best done when the source is highly mag- 
nified as it crosses a caustic, as this minimizes the effects of 
blended light and increases the signal-to-noise. Thus such 
measurements require target-of-opportunity observations. 
For highly-magnified events, spectra with S/N <~ 100 per 

9 Note that blending in generally not a problem, as binary-lens events can typically be easily deblended. 

10 For online alerts, see http://www.astrouw.edu.pl/~ftp/ogle/ogle3/ews/ews.html (OGLE) and http://www.roe.ac.uk/~iab/alert/alert.html 
(MOA). 



FIG. 3.- The cumulative distribution of I magnitudes for OGLE and 
MOA alerts in 2002. The solid line is for OGLE, the dotted line is 
for MOA, and the dashed line is the total. The majority (~ 90%) of 
alerts are from OGLE. 

spectral typing, high resolution spectra are desirable for a 
number of other applications, including resolution of the 
atmosphere of the source star (Gaudi & Gould 1999, Cas- 
tro et al. 2001, Albrow et al. 2001b), detailed abundance 
analysis (Minniti et al. 1998), and detection of a luminous 
lens (Mao, Reetz, & Lennon 1998). Note also that as a 
byproduct, true space velocities of a sample of stars in the 
bulge will be obtained by combining the proper motions 
and parallaxes of the sources acquired from astrometric 
measurements with radial velocities determined from the 
spectra. 

4.5. An Example Campaign 

In this section, we review the requirements for measur- 
ing 6* for the source stars of Galactic bulge microlensing 
events, and outline the resources needed for a campaign 
aimed at measuring for a significant number of sources. 

The first requirement is a large sample of caustic- 
crossing binary-lens events from which to choose targets, 
which in turn requires an even larger sample of microlens- 
ing events. A large sample is important in that it ensures 
that only interesting and promising sources and events are 
followed. Currently, both the OGLE and MOA collab- 
orations monitor many millions of stars in the Galactic 
bulge. Both reduce their data real-time, enabling them to 
issue 'alerts,' notification of ongoing microlensing events 
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(Udalski et al. 1994, Bond et al. 2002a). 10 Combined, 
these two collaborations should alert about 500 events per 
year (with the majority of alerts from OGLE). Extrapo- 
lating from previous results (Alcock et al. 2000a, Udalski 
et al. 2000), approximately 5% of these will be caustic- 
crossing binaries, or 25 events per year. Figure 3 shows 
the cumulative distribution of apparent (i.e. uncorrected 
for rcdcnning) I magnitudes of the 438 independent bulge 
microlensing alerts in 2002 for which baseline magnitudes 
were available. Of these events, 382 (87%) were alerted 
by OGLE, 61 (14%) were alerted by MOA, and 5 were 
alerted by both collaborations. For the typical colors of 
sources in the bulge (V — I ~ 2), only about 50% of the 
alerts have V < 20, and thus would have been accessible 
to SIM. The boundary between dwarfs and giants will oc- 
cur at an apparent magnitude that depends on the color of 
the source, the distance to the source, and the reddening. 
However, for definiteness we will simply assume that the 
boundary between giants and main sequence stars occurs 
at roughly / = 17. With this assumption, we can there- 
fore expect that approximately 20% of all alerted events 
will be due to giant stars. Therefore, we can expect ap- 
proximately 20% x 5% ~ 1% of all alerts, or ~ 5 events 
(assuming 500 alerts), to be caustic crossing events with 
giant sources, and 80% x 5% <~ 4%, or <~ 20 events, to be 
caustic crossing events with main sequence sources, ~ 8 of 
which will be bright enough to monitor with SIM. 11 These 
numbers are likely to remain valid at least for the next 
several years. However, in the more distant future, and 
in particular by the time SIM is launched, it is likely that 
the next generation of microlensing survey collaborations 
will have come online. Thus we can expect that, when 
SIM time is operational, 12 a considerably larger sample of 
caustic-crossing events will be available. 

Survey-type experiments are needed to discover mi- 
crolensing events toward the bulge, and survey-quality 
data is generally sufficient to uncover the caustic-crossing 
nature of the target events. However, as we discussed in 
§3.1, more accurate and densely-sampled photometry is 
generally needed during the caustic crossing in order to 
measure At. Currently, there are several collaborations 
with dedicated (or substantial) access to l-2m class tele- 
scopes distributed throughout the southern hemisphere 
that closely monitor alerted microlensing events with the 
goal of discovering deviations from the single-lens form, 
with emphasis on search for extrasolar planets (Albrow 
et al. 1998, Rhie et al. 2000, Tsapras et al. 2001, Bond 
et al. 2002b). These collaborations have also been quite 
successful in predicting and monitoring binary-lens caus- 
tic crossings. It seems likely that these collaborations, or 
similar ones, will still be in place when the next generation 
of interferometers, or even SIM, come online. 

In our Monte Carlo simulations we derived the expected 
precisions <rg E /0 E assuming that the photometric errors 
were dominated by photon statistics, and that a total of 
N = 60, 000 photons where collected over the entire ex- 



posure time for each event. This corresponds to total 
exposure time of T = 1.6 hour for SIM on an J = 18 
source. We assume the fractional photometric uncertainty 
is A -1 / 2 , and that the astrometric uncertainty was re- 
lated to the photometric uncertainty via the expression, 
erg = N~ 1 > 2 8f, with Of — 2.5mas, as appropriate for SIM. 
Since we assumed that the photometric and astrometric 
errors are given simply by photon statistics, it is trivial 
to scale our results for other total exposure times T and 
source brightnesses assuming the characteristics of SIM: 
<T0 E /6 E oc T- 1 / 2 , and ae E /9 E oc io°- 2 ( J - 18 ). For the pur- 
poses of planning observations, and providing an order-of- 
magnitude estimate for the number of angular radii that 
can be measured for a given amount of SIM time, it is 
useful to derive an expression for the exposure time re- 
quired to achieve a given median photometric precision. 
To be conservative, we assume that the blending is small, 
but non-negligible. Specifically, we adopt the median error 
found for the Monte Carlo simulations assuming a blend 
fraction of f t = 1%, which is ag E /8 E ~ 5.5%. Then, 



2 hr 



5% 



10 0.4(/-18) 



(15) 



Thus, for giant sources (/ ~ 15), 7.4 minutes are required 
to achieve 5% precision, whereas for main-sequence sources 
(7 ~ 20), 12.3 hours are required for 5% precision, or 1.4 
hours for 15% precision. 

We have focussed here primarily on astrometric obser- 
vations with SIM because its capabilities are well-suited 
to this application. The basic requirements to be able to 
measure E accurately for the events we have discussed are 
relatively high astrometric precisions, ~ 10 fias, and high 
sensitivity (via, i.e. large apertures), as the sources we are 
considering are faint, / = 15 — 20. These faint sources 
are inaccessible to current ground-based interferometers. 
Upcoming large-aperture, ground-based interferometers, 
such as the Very Large Telescope Interferometer or the 
Keck Interferometer, should be able to achieve the requi- 
site astrometric precisions on all of the bright (I <J 15) 
giant events. If the target microlensing source happens 
to have a bright star within the isoplanatic angle, it may 
be possible to employ phase referencing to extend sensi- 
tivity to very faint (I <; 20) sources. This would allow 
one to measure E for main-sequences sources from the 
ground as well. Finally, it may be possible to determine 
#e from single-epoch measurements of the visibility and/or 
closure phase (Delplancke, Gorski, & Richichi 2001, Dalai 
& Lane 2003). In this way, sensitivity could plausibly be 
extended to main-sequences sources by making carefully- 
timed interfcrometric measurements of the source when it 
is highly-magnified during a caustic crossing. However, it 
is not clear if there exists enough structure in the image po- 
sitions during this time to extract E . This remains an in- 
teresting topic for future study. Non-targeted space-based 
astrometric surveys, such as the Global Astrometric Inter- 
ferometer for Astrophysics, are generally not well-suited 

11 This assumes that the detection efficiency of binary-lens events does not depend on the /-magnitude of the source. In fact, deviations from 
the single-lens form will generally be easier to detect in brighter sources, however this bias is likely to be relatively small for caustic-crossing 
events, which generally exhibit dramatic and easily-detectable deviations from the point lens form. This may seem in contradiction with the 
fact that none of the five events toward the bulge presented in Table 1 are main sequence sources. However, this almost certainly a selection 
effect: bright binary-lens events are currently preferentially monitored by the follow-up collaborations, in order to achieve higher signal-to- noise 
during the second caustic crossing. 

12 The projected launch date for SIM is currently 2009. 
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to this application, due to the relatively sparse sampling 
of the source stars. 

Finally, access to target-of-opportunity time on 8-10m 
class telescopes would allow for accurate spectral typing 
of the source stars. Several nights per bulge season would 
likely be adequate to type the <~ 25 caustic-crossing events 
per year. However, more time would be required to per- 
form some of the auxiliary science discussed in §4.3, such 
as resolution of the source-star atmospheres. 

Thus, by combining alerts from survey collabora- 
tions, with comprehensive ground-based photometry from 
follow-up collaborations with access to dedicated (or semi- 
dedicated) lm-class telescopes, and a modest allocation of 
a total 10 hours of SIM time, it should be possible to mea- 
sure the angular radii of ~ 80 giant stars in the bulge to 
5%, or ~ 7 main sequence stars to 15%. Several nights of 
target-of-opportunity time on 8 — 10m telescopes should 
allow for accurate spectral typing of the sources via high 
or low-resolution spectroscopy. 

5. CONCLUSION 

We have outlined a method to measure the angular radii 
6* of giant and main sequence source stars of fold caustic- 
crossing binary microlensing events toward the Galactic 
bulge. Our method to measure 6* consists of four steps. 
First, survey-quality data can be used to discover and alert 
caustic-crossing binary-lensing events. Such data is suffi- 
cient to characterize the event timescale t& and the angle 
4> of source trajectory with respect to the caustic. Dense 
sampling of one of the caustic crossings yields the caustic- 
crossing timescale At. The global solution to the binary- 
lens light curve yields a prediction for the trajectory of the 



centroid of the source up to an unknown angle a, and the 
scale, 6e- Thus a few, precise astromctric measurements 
during the course of the event yield 6*e- The angular source 
radius is then simply given by 0* = sin0. 

We argued, based on past experience with modeling 
binary-lens events, that the parameters At, 0, and £e 
should be measurable to a few percent accuracy, provided 
one caustic-crossing is densely and accurately sampled, 
and the entire event is reasonably well-covered. 

We then performed a series of Monte Carlo experiments 
that demonstrated that astrometric measurements during 
the course of the binary-lens event should allow for the 
determination of 6*e to ~ 2% accuracy, assuming photon- 
limited statistics and a total of 60,000 photons per event. 
This is a factor of ~ 50 fewer photons than are required 
to measure 9e to the same precision in single-lens events 
and corresponds to an exposure time of T = 1.6 hour with 
SIM on an I = 18 source. Therefore, it should be possible 
to measure 6* for a significant sample of giant and main- 
sequence stars in the bulge with reasonable expenditure of 
resources. 



We would like to thank Neal Dalai for helpful conversa- 
tions. We would also like to thank the anonymous referee 
for useful comments and suggestions. This work was sup- 
ported by NASA through a Hubble Fellowship grant from 
the Space Telescope Science Institute, which is operated by 
the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555, by JPL contract 
1226901, and by the Science Research Center (SRC) of the 
Korean Science and Engineering Foundation (KOSEF). 



REFERENCES 



Afonso, C. ct al. 2000, ApJ, 532, 340 
Afonso, C. et al. 2001, A&A, 378, 1014 
Albrow, M. et al. 1998, ApJ, 509, 687 
Albrow, M. D. et al. 1999a, ApJ, 522, 1011 
Albrow, M., et al. 1999b, ApJ, 522, 1022 
Albrow, M. D. ct al. 2000, ApJ, 534, 894 
Albrow, M. D. et al. 2001a, ApJ, 549, 759 
Albrow, M. ct al. 2001b, ApJ, 550, L173 
Albrow, M. D. ct al. 2002, ApJ, 572, 1031 
Alcock, C. ct al. 1997, ApJ, 491, 436 
Alcock, C. ct al. 2000a, ApJ, 541, 270 
Alcock, C. ct al. 2000b, ApJ, 542, 281 
An, J. A. et al. 2002, ApJ, 572, 521 
Armstrong, J. T. et al. 1998, ApJ, 496, 550 
Baltz, E. A. & Gondolo, P. 2001, ApJ, 559, 41 
Bond, I. A. et al. 2001, MNRAS, 327, 868 
Bond, I. A. et al. 2002a, MNRAS, 331, L19 
Bond, I. A. et al. 2002b, MNRAS, 333, 71 

Castro, S., Pogge, R. W., Rich, R. M., DePoy, D. L., & Gould, A. 

2001, ApJ, 548, L197 
Chang, K. & Refsdal, S. 1979, Nature, 282, 561 
Chang, K. & Refsdal, S. 1984, A&A, 132, 168 
Colavita, M. M. et al. 1999, ApJ, 510, 505 
Dalai, N. & Griest, K. 2001, ApJ, 561, 481 

Dalai, N. & Lane, B. F. 2002, ApJ, submitted (astro-ph/0209249) 

Dclplancke, F., Gorski, K. M., & Richichi, A. 2001, A&A, 375, 701 

Dominik, M. 1998, A&A, 329, 361 

Dominik, M. 1999a, A&A, 341, 943 

Dominik, M. 1999b, A&A, 349, 108 

Gaudi, B. S. & Gould, A. 1999, ApJ, 513, 619 

Gaudi, B.S., & Petters, A.O. 2002a, ApJ, 574, 970 

Gaudi, B.S., & Petters, A.O. 2002b, ApJ, 580, 468 

Gould, A. 1994, ApJ, 421, L71 

Gould, A. 2000, ApJ, 535, 928 

Gould, A. 2001, PASP, 113, 903 

Gould, A., & Salim, S. 1999, ApJ, 524, 794 

Gould, A., & Han, C. 2000, ApJ, 538, 653 

Gould, A., & Welch, D. L. 1996, ApJ, 464, 212 



Graff, D.S., & Gould, A. 2002, ApJ, 253 
Han, C. & Kim, T. 1999, MNRAS, 305, 795 
Han, C, & Kim, H. 2000, ApJ, 528, 687 
Han, C, Chun, M., & Chang, K. 1999, ApJ, 526, 405 
Jaroszyhski, M., & Mao, S. 2001, MNRAS, 325, 1546 
Jeong, Y., Han, C, Park, S.-H. 1999, ApJ, 511, 569 
Lcnnon, D. J., Mao, S., Fuhrmann, K., & Gehren, T. 1996, ApJ, 471, 
L23 

Mao, S. & Witt, H. J. 1998, MNRAS, 300, 1041 

Mao, S., Reetz, J., & Lennon, D. J. 1998, A&A, 338, 56 

Minniti, D., Vandehei, T., Cook, K. H., Griest, K., & Alcock, C. 

1998, ApJ, 499, L175 
Nordgren, T. E. et al. 1999, AJ, 118, 3032 
Paczyhski, B. 1998, ApJ, 494, L23 

Petters, A. O., Levine, H., & Wambsganss, J. 2001, Singularity 

Theory and Gravitational Lensing (Boston: Birkhauser). 
Popper, D. M. 1998, PASP, 110, 919 

Rhie, S.H. & Bennett, D. 1999, preprint (astro-ph/9912050) 
Rhie, S. H. et al. 2000, ApJ, 533, 378 

Richichi, A., Fabbroni, L., Ragland, S., & Scholz, M. 1999, A&A, 
344, 511 

Sahu, K. C. 1994, Nature, 370, 275 

Schneider, P., Ehlers, J., & Falco, E. E. 1992, Gravitational Lenses 

(Berlin: Springer) 
Schneider, P. & Weiss, A. 1986, A&A, 164, 237 
Schneider, P. & Weiss, A. 1987, A&A, 171, 49 
Schneider, P. & Weiss, A. 1992, A&A, 260, 1 
Tsapras, Y. et al. 2001, MNRAS, 325, 1205 

Udalski, A., Szymanski, M., Kaluzny, J., Kubiak, M., Mateo, M., 
Krzeminski, W, & Paczynski, B. 1994, Acta Astronomica, 44, 227 

Udalski, A., Zcbrun, K., Szymanski, M., Kubiak, M., Pietrzynski, 
G., Soszynski, I., & Wozniak, P. 2000, Acta Astronomica, 50, 1 

van Belle, G. T. 1999, PASP, 111, 1515 

van Belle, G. T., et al. 1999a, AJ, 117, 521 

Wozniak, P. & Paczyhski, B. 1997, ApJ, 487, 55 

Wozniak, P. R., Udalski, A., Szymanski, M., Kubiak, M., Pietrzynski, 
G., Soszynski, I., & Zebrun, K. 2001, Acta Astronomica, 51, 175 



